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Abstract. 

We propose a generalization of the replica trick that allows to calculate the 
large distance asymptotic of off-diagonal correlation functions in anyonic models 
with a proper factorizablo ground-state wave-function. We apply this new method 
to the exact determination of all the harmonic terms of the correlations of a gas 
of impenetrable anyons and to the Calogero Sutherland model. Our findings are 
checked against available analytic and numerical results. 
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1. Introduction 

Generalized anyonic statistics, which interpolate continously between bosons and 
fermions, are considered one of the most remarkable breakthrough of modern physics. 
In fact, while in three dimensions particles can be only bosons or fermions, in lower 
dimensionality they can experience exchange properties intermediate between the two 
standard ones [1]. In two spatial dimensions, it is well known that fractional braiding 
statistics describe the elementary excitations in quantum Hall effect, motivating a 
large effort towards their complete understanding. 

In one dimension, anyonic statistics are described in terms of fields that at 
different points {xi ^ x^) satisfy the commutation relations 

^\{x^)^\{x^) = e--^(-i--^)4'^(a;2)*^(,ri) , (1) 

'^a{xi)^\{x2) - e-^"'^<^'-^-^^\{x2)^A{xi) , 

where e{z) = —e{—z) = 1 for 2; > and e(0) = 0. k is called statistical parameter and 
equals for bosons and 1 for fermions. Other values of k give rise to general anyonic 
statistics "interpolating" between the two familiar ones. 

After the second life of one-dimensional gases started following their experimental 
realization in cold atomic setups [2], there have been several proposal to create and 
manipulate anyons [3], even in one dimension. A few ID anyonic models have been 
introduced and investigated [4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 
22, 23, 24, 25, 26]. In this paper we consider the calculation of off-diagonal correlation 
(or one-particle density matrix) for a system of N particle 

gnx)=p%{x) = {'^':^ix)nm, m 

of models that have a factorizable many-body ground-state wave function as 

...,xn) = Cn'[[ /(e^^"''/^ - e2--'/^)B,(a;fc - xi) , (3) 

k<l 

with f{x) a generic even function. Cat a normalization constant, and 

f e""'' a; < 0, 
B4x) - r ^ ; 4 

II X > 0. 

that completely specify the statistic. Among these models, the more relevant are the 
anyonic generalization of Tonks-Girardeau gas and of the Calogero-Suthcrland model. 

In this paper we propose a generalization of the replica trick that is able to select 
the correct anyonic branch in the replicated space. The use of this replica method 
started in the context of the diagonal correlation of the Calogero-Sutherland model 
[27], has been later extended to a gas of impenetrable bosons [28, 29], and finally to 
the off-diagonal correlation of Calogero-Sutherland, but only for bosonic and fermionic 
statistics [30]. It has been already proposed [28] that in the replica approach one 
deals with two functions of the replica number n, one for fermions and one for bosons, 
or better a single function with two branches, a fermionic and a bosonic ones that 
coincide only for integer n. Here we push forward this interpretation and we claim 
that the replicated correlation has an infinite number of branches corresponding to 
the different anyonic statistics. The apparent lack of rigor of this approach (that is 
shared by most of replicas calculations starting from the celebrated Parisi mean-field 
solution of the spin-glass [31]) can be justified a posteriori by the amazingly simple 
final results that agree with all the available analytics and with the direct numerical 
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calculations. The analytic continuation we propose in this paper allows us to obtain 
exact asymptotic expansions for large distances that are valid for the off-diagonal 
correlation for general anyonic parameter k. We also exploit the property that the 
anyonic correlation function satisfies a Painleve V differential equation [21] , that allows 
a rigorous analytic check of the results obtained trough replicas. To further check the 
replica predictions, we perform a careful numerical comparison of the finite-size off- 
diagonal correlation obtained via the Toeplitz determinant representation [8, 21]. 

The manuscript is organized as follows. In Sec. 2 we introduce the anyonic Tonks- 
Girardeau gas and the strategy for the replica calculation. In Sec. 3, starting from the 
example of a simple integral we show how to select the correct anyonic statistic from 
the replicas. In Sec. 4 we review the replica calculation for the bosonic impenetrable 
gas following Ref. [28]. In Sec. 5 we apply the replicas to the anyonic Tonks- Girardeau 
gas, showing that all the harmonic terms in the large distance expansion are captured 
by the saddle-point approximation in the replica space. In Sec. 6 we calculate the 
corrections to the saddle-point approximation, both by perturbation theory and via 
the solution of the Painleve V equation. In Sec. 7 we compute numerically the off- 
diagonal correlation function, showing perfect agreement with the analytic results. In 
Sec. 8 we use the replica trick to obtain the correlation function in the fully anyonic 
Calogero-Sutherland model. Finally in Sec. 9 we draw our conclusions. 



2. The anyonic Tonks-Girardeau gas 



In first-quantization language, the Lieb-Liniger Hamiltonian is 

i * l<i<j<N 

where the 7V-anyons ground-state function exhibits the generalized symmetry under 
the exchange of particles. For k = the model reduces to the bosonic Lieb-Liniger [32], 
while for k = 1 to free fermions. The Tonks-Girardeau gas is obtained for c oo and 
corresponds to impenetrable anyons. In this limit, the ground-state wave function has 
the form of Eq. (3) with f(y) = \y\. This anyonic model has been explicitly studied 
in Refs. [4, 7, -5, 8, 11, 12, 22, 21, 24]. 

The anyonic one-particle density matrix for -|- 1 particles can be written [8], as 
A''-dimensional integral (in the variables Xj/L = tj and t = x/L) 

TV 

Y[B^{t,-t)\sm[TT{ts-t)]\\sm[Trts]\ 4sin[7r(tj - fi)]". (6) 

s=l l<i<j<N 

Using the identity 

n 4sin[7rfe-tfe)]'= n |e^'^*'■-e^'-*^|2 = |A^({^0)|^ (7) 

\<j<k<N \<j<k<N 

we can identify the second product in the integral (6) with the square of the absolute 
value of a Vandermonde determinant. Thus the reduced density matrix is an average 
over the distribution [AtvI^: 

P%+r{t) = {\{B^{ts-t)\z-Zs\\l-Zs\\ , (8) 

\«=1 / |A„|2 
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where Zg = e'^'^*" and z = e'^'^*. 

In the case of bosons (k = 0) the reduced density matrix can be obtained [28] 
through a replica trick considering for a positive integer n 

ZS='{t) = (f[{z-z,r"{l-z,)A , (9) 

and then p'lfj^i{t) is obtained from Z2n{t) by a suitable analytic continuation to 
n — > 1/2. It has been shown by Kurchan [33] that to obtain the average of the 
absolute value from the replicated correlation, we have to require that the analytic 
continuation for any complex n with positive real part should diverge as maximum as 
a power law. 

How to modify such replica calculation to anyons when we deal with averages 
over a function of unitary modulus but with arbitrary phase? Our trick is to start by 
considering an anyonic parameter that is a rational number, i.e. k = q/p with p and 
q integer. We can then consider 



= 1 



|AnP 



that, after a proper analytic continuation, for 2np 1 could give the desired reduced 
density matrix. The ambiguity in this case is obvious, in fact Z2np includes all the 
correlations of anyons with statistical parameter k = q/p for any q = 1 . . . 2p. This 
means that different analytic continuations will give different anyonic correlations 
including bosonic and fermionic ones. We will propose a general criterion that would 
allow to choose the proper analytic continuation for any value of the statistical 
parameter. To elucidate the meaning of our trick we start by applying it to a trivial 
integral that contains all the essence of the modified replica trick. 

3. The replica trick for an integral 

Let us consider the integral 

/ = / dt\t\e-* =1, (11) 



that has the replica representation 

/+0O 1 I 2nin 
dte^e-' r(l/2 + n). (12) 
-oo ^ 

Notice that Ii = 0, while considering n integer and analytically continuing we have 
J|„ = r(l/2 + n), that leads to /|„ = 1, the desired result. This analytic continuation 
selected the "bosonic" branch of the integral, giving the integral with the absolute 
value, because it is the only analytic continuation that does not diverge exponentially 
in the complex plane for n ioo, i.e. that satisfy Kurchan criterion for the absolute 
value [33] . Oppositely taking directly I^n we would have obtained the "fermionic" 
branch. 

Let us then modify the integral anyonically: 

r+oo , , , , _.2 1 + e-^'^'/P 



r 



dtB,{t)\t\e-' =^^^ , (13) 

-oo ^ 
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with Bi^{t) given by Eq. (4) and k = q/p. The proposed analytic continuation is 

dtt^^Pe-' T{l/2 + np), (14) 

-oo ^ 

where wc loose any knowledge of the value of q, that is the previous mentioned 
ambiguity. Considering n and p integers and analytically continuing one would always 
get the bosonic branch. On the other hand, one can play the following game 

1 I „-2iTTnq„2iTm{p+q) 

hup = — ^ r(l/2 + np) , (15) 

with q, p, n integers. Considering the analytic continuation 

1 _|_ p-2i-Knq 

lLp = —-^ r(l/2 + np), (16) 

we get 

1 I p-i-Kq/p 

n = , (17) 

that selects the desired anyonic branch. The reason why this simple game works is 
that /jnp is the only analytic continuation which has the correct asymptotic behavior 
for large imaginary n (and positive real part) that selects the correct branch. The 
message is that the integral /2„p contains all statistics from bosons to fermions and a 
given one can be obtained just by selecting the correct behavior at large imaginary n. 

With this observation in mind, we can then use the result Z2n (t) for the bosonic 
Tonks-Girardeau gas and for the Calogero-Sutherland model. Properly continuing 
Z2np{t), we easily get the anyonic result. Thus in the next section we review the 
bosonic result and in the following one we generalize it to anyons. 

4. Bosons review 

All the content of this section is taken from the pioneering work of Gangardt [28] . We 
report this material here only to make this paper self-contained. 

The replicated average for bosons Eq. (9) is (we recall z = e^'^**) 

Z2n{t) = f'^ d^t|A;v(e''^^*)|' n 11 + e'"'*l'"|^ + e^-*^!^", (18) 

where we have normalized .^2n(l) = 1) introducing the constant MN{2n) (see Ref. 
[28] for its precise value, it is not essential for what follows). The integral (18) has a 
dual representation [34] by an integral over n variables 

Nn rl 2n 

Z2n{t) = ^ / d'-X Al^ix) n (1 - (1 - ^)^a)^, (19) 
S2n Jo 



where the normalization constant is given by the Selberg integral 

r2(a)r(l + a) 



'2n 



.1 2n 

/ d'-xAUx)^l[ 

Jo „_i 



, r(2n + a) 

a— 1 ^ ' 



(20) 



In Eq. (19) the number of particles N appears only as a parameter. This 
representation allows us to obtain the asymptotic expression for Z2n suitable for 
analytic continuation in n. In the large A'' limit the integrand in (19) oscillates rapidly 
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and the main contribution conies from the endpoints x± = 1,0 which are the only 
stationary points of the phase. We change variables near each endpoint 

= T77T^ — ^' a=l,...,l, 

jV(l-^-i) ' b = l + l,...,2n, (21) 

The integrand in (19) simplifies in the large N limit: 



and the integration measure including the Vandermondc determinant is factorizcd as 

d'aAf(gJd^"-'6AL-,(a) 
[7V(l-^)p[7V(l-0-i)](2"-O= 

The remaining integrals are calculated using [36, 37, 38] 



^00 ^ J- 

i{x)= d'^aAuu n = ^"'^ n 

"'O a=l a=l 



+ a). (24) 



Multiplying the contribution of each saddle point by the number of ways to distribute 
variables we obtain the asymptotic of the integral (19) as a sum of 2n + 1 terms | 



where we have introduced X = N{z^/^ — z^^/^)/2i = TVsinTrt and the factors 

pi -n _ G(l + l)G(2n-^ + l) 

^2„-ll p(2n-Z + a) " G(2n + 1) ' ^ 

where G{x) is the Barnes G-function 

oo 

G{z + 1) = (2,r)^/2e-(^+(^-+i)^')/2 JJ (1 + z/k)'e-'+'"/^'"'^ , (27) 

fc=i 

7b is the Euler constant, G{1) = G(2) = 1, and G(3/2) = 1.06922.... Using the 
Barnes fimction is particularly useful to get the analytic continuation easily. Notice 
that also at this level, if we set 2n = 1 in Eq. (25) we still get the fermion result 
3«=i(a:) = sin((iV + l)Trt)/{{N + 1) siuTrt). 

For n integer, we can extend the sum from < ^ < 2n to all the integers, because 
vanishes for i < and I > 2n (it is a simple properties of the Barnes G function 
inherited from the F). While this is an innocuous change for integer n, it is the 
fundamental step that will allow for replica symmetry breaking, because for general 
values of n (i.e. non integers) are non zero for any I and the sum results in a true 
infinite series. Changing the summation variable to fc = Z — n we obtain finally 

^^"(*) - {2N\sinnt\rr.^ ^^'^^ + ^ ^(-1) [F,^ \ ^^j^, ^^2 j 



I we corrected an error/typo irrelevant in the limit n — ► 1/2 in Ref. [28], i.e. a factor that is one 
when n — > 1/2. 
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This apparently useless shift is the crucial point where the replica symmetry is finally 

broken because it selects the zero mode. Now we arc in a position to take the limit 
n ^ 1/2 which results in the desired harmonic expansion of the one-body density 
matrix 



PN+lit) 



1 27V sin TTt 1 1/2 



f:(-i)'iG(3/2+.)G(3/2-m' -^;y;^y 



The form of this expansion has been firstly predicted by Haldane [35] , but the complete 
knowledge of all the amplitude has been only obtained by means of replicas [28] . 



5. The anyons 

In the calculation of Z2np everything proceeds like in the calculation for boson, leading 
to the replacement n ^ np in the replicated correlation (25) 



^ rClrrn + ^ \ F.I f AN+2np){np-l) 



(29) 



The correct analytic continuation to get the anyonic branch with k = q/p is obtained 
by selecting the behavior at large imaginary n. As for bosons, this is obtained by 
choosing the zero mode in the sum over I and then by extending the sum over all 
integers. Thus we set 

np — I = nq + m, (30) 

that, amazingly, is the only change we need to describe anyons from the bosonic 
formulae. With this substitution we get 



(31) 



Y(c) ^ C2X)("P^("9+'"))^ + ("P+("9+'"))2 

that is ready for the analytic continuation 2np 1, obtaining 

oo '_pV2-9/2p-m ^ ^(Ar+l)(g/2p+m) 

Pn+A^) = X] (-l)^^/^^"*"™' (■2XlV2+2(r»+<//2p)2 ' 

(32) 

where only the ratio k = q/p appear and so it can be generalized to any k not only 
ratio of integers, giving (restoring z = e^^trx/L j 

'3 _ K _ rn)n(^ -I- £ -I- m)] ^ f,i-iTx(K+2m){N+l)/L 

.(33) 



o« (X)- V f iy./2+^) [G(|-t-m)G(| + f+m)]^ e- 



(2iV| sin7rx/L|)i/2+2('"+«/2)^ 

According to bosonization (or harmonic fluid approach), any anyonic off-diagonal 
correlation in the thermodynamic limit N ^ oo with po = N/L kept constant, admits 
the expansion [11] 

giix)=popoo 5Z (2^+^)2 —■ (34) 

Here d{x) = |Lsin(7rx/I/)|, are unknown non-universal amplitudes, while K is 
the universal Luttinger liquid exponent equal to 1 for Tonks-Girardeau (we fixed 
bo = 1). The two expansions agree perfectly, confirming the correctness of the result. 



Off-diagonal correlations in ID anyonic models 



8 



Furthermore, we obtained an exact expression for all the coefficients in the harmonic 
expansion 



Poo ^77 



[G(|-f-m)G( 



f +m)]' 



with Poo = 



[G(|-f)G(| + f)]^ 



(35) 



2l/2+2(m+K/2)2 ..^-"f'oo ^ 

The leading term poo has been already obtained from Fisher-Hartwig conjecture 
applied to the Toeplitz determinant [8, 21], and the two results agree. All the other 
terms are new. It is important to write down the ratio of two consecutive amplitudes 

T(3/2 + ^/2 + m) 1%i+2.+2m .o.^ 

r(l/2-«;/2-m)J ' ^ ' 

The first corrections to the leading behavior are given by the terms with m = ±1, 
that satisfy 

M-. = 4(l-.^r. (37) 

We will see that the same result also follows from the Painleve V equation, given a 
further check of the replica approach. 



6. Beyond the saddle points 

The formula just derived gives the leading term in 1/N of the desired correlation 
function. As well known, this does not give all the dominating terms in the 
thermodynamic limit, defined as N ^ oo with Nnx/L = kpx taken constant. One 
then needs to calculate corrections to the previous terms. A systematic way to get 
these corrections is to perform the perturbation theory close to the saddle points order- 
by-order considering in Eq. (21) the full expansion in ^a,b close to x± = 0, 1, as done 
in Rcf. [28] for bosons. This approach is very general and straightforward and gives 
all the 1/N corrections, but unfortunately becomes soon computationally cumbersome 
increasing the order. For this reason, we only report the first-order computation in the 
next subsection. Alternatively, we can take a different approach based on the solution 
of the Painleve V equation, that gives very easily all the terms in the expansion, but 
only in the thermodynamic limit. 



6.1. First order in perturbation theory 

In Ref. [28] the first order correction in to Z2n{t) has been explictly calculated (see 
Eq. (65) there). Using this result and simply posing 2n 2np and then performing 
the shift np — I = nq + m, we get for the first order term (we recall that X = A'' sin irt 

so it is 0{N)) 

— {np — nq — m){np -\-nq-\- m)[{np + ng + m)e^'^^ — {np — nq — m)e~^'^^] 

- [{np + nq + mf + {np-nq- mf] , (38) 

taking 2np — > 1 and then considering only the zero mode m = (higher modes give 
rise to higher corrections in 1/N), we get 

1 1 - 1 

-7777(1 + '«^)+« K— T-r. (39) 

2N^ ^ 4 Nt&niTx/L ^ ' 




Notice that for bosons and fermions {k = 0, 1) the tan correction is absent: the role 
of this term in the rephca approach is only understood while studying anyons. The 
first term — (1 + n^)/2N just corresponds to the natural replacement — > + 1 in 
the denominator of the leading term. 



6.2. Large distance asymptotic expansion via Painleve V equation 

In Ref. [21] we have shown that pK.{t) for finite N satisfies second order differential 
equation of Painleve VI type. The differential equation does not depend on the 
anyonic parameter, i.e. it is the same for any k. What characterizes the statistic 
are the boundary conditions at t = that do depend on the anyonic parameter. It is 
well known that in the thermodynamic limit, the Painleve VI equation reduces to a 
Painleve V (see e.g. [36, 39, 40])) that has been firstly derived for bosons by Jimbo et 
al. [41] 

ixa"f + A{x(7' - 1 - CT){xa' + cr'^ - cr) = , (40) 

•f^ t \ d\ogp{x) 

with <y[x) — X ; . (41) 

dx 

To get a systematic expansion for large x from this equation we need to impose the 
anyonic boundary conditions at x = oo. We need the knowledge of the first two leading 
terms in the harmonic-fluid expansion 6o ^^nd bi. Actually 6o gives only the overall 
normalization, in fact (j{x) does not depend on it, while 5_i is essential to fix all the 
signs (in fact an old sign error in 6_i = 6i for bosons [42] caused all the expansion in 
[41] to have the wrong signs. This has been corrected only in [28], see also [43]). 



Off-diagonal correlations in ID anyonic models 



10 



In this subsection we set kp = 1 for simplicity in writing, the correct formulas 

are obtained by replacing x kpx. The complete asymptotic expansion of p{x) 
is obtained by multiplying any term in the harmonic expansion (34) by an analytic 
function at a; = oo, i.e. 



p{x) 



X 



1/2+K2/2 



1 + — + 4 



\ X X'^ J 



, e ™ ci C2 \ 



(42) 



Form this we have the large distance expansion of a that fixes the boundary conditions: 

lim a{x) = inx - - 2i6_ie-2".T-i+2'^ + 0{x°) , (43) 

where o is some not better specified exponent that will follow from the equation. 
Plugging this expansion in the Painleve equation one can get iteratively all the various 
coefficient Oi, 6;, c;. rf;. The first term in this expansion is Eq. (37), that provides a 
further consistency check of the replica approach. After long and tedious algebra we 
got the first coefficients: 

1 _ ^2 

rfi = , (44) 

c,= -iii±^Hi±^fc^, (45) 
a,=z(i^*^^fc^, (46) 
rf2 = ^^(«^ + 4«2-l), (47) 
02 = - ^(3 - k){1 - k)(31 - 48k + 28^^ - 8k^ + k^) , (48) 

C2 = - ;^(3 + k)(1 + k)(31 + 48k + 28k2 + 8/«-^ + k*), (49) 

that for K ^ reproduce the known result for bosons [28, 41]. Notice that d\ 
reproduces the result from first order perturbation theory. In Fig. 1 we report all 
these coefficients as function of k. Notice that since at the fermionic point all the 
terms but two must vanish we have = dm = ^ for any m, while the vanishing of 
the other terms is ensured by 6i = 0. At the bosonic point instead the reality of the 
correlator imply ai = c* that is satisfied by our expressions. 

7. The comparison with numerics 

In this section we show how our expansions perfectly fit with the numerical data 
obtained by the exact evaluation of the off'-diagonal correlators at finite N = L through 
the evaluation of the Toeplitz determinants derived in Refs. [8, 21]. 
At finite volume and unit density, the expansion can be written as 

p{x) = 



1[1 + o{N-')] + ^-^hzTT + ..rJL.ru + 0{N-') 



Ntamrx/L [A^sin7r.7;/_L 



[A^sinTTx/L] 2+ 2 

^-2iNx f ai a2 3 ^ 

^''"^[Arsin7ra;/L]2-2« y-'^ NUnnx/L ^ [Nsimvx/L]-^ ^ \ 
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Figure 2. Off-diagonal correlation function for k = 0.1 and N = 81. Top panels: 
Real (left) and imaginary (right) parts of Pgi(t). Second rows: p^'^ {t)/ ppH{t) — l. 
Third rows: Subtraction of the first two harmonic terms b±i. Fourth rows: 
Further subtraction of the analytic contributions di 2- Bottomost panels: Further 
subtraction of the analytic corrections to the harmonic terms ai and ci. 



-hi. 



[N sin7rx/L]2+2'^ 
+higher harmonics 



Cl 



C2 



N tamrx/ L [Nsimrx/L]-' 



0{N 



-3< 



(50) 



Notice that we have used in the leading term N instead of iV — 1 to cancel the 1 /N 
correction to the constant in the first line, according to Eq. (39). It has been shown 
already in Ref. [8] that only the leading term in this expansion (obtained from Fisher- 
Hartwig conjecture) almost perfectly describe p{x) for all k < 1/2. In Ref. [11] it has 
been argued that for larger value of k the first harmonic term (i.e. 6_i) is fundamental 
to have a good description of the asymptotic results because its power 2 — 2k gets closer 
to zero. Including this term, we have a perfect description of the asymptotic behavior 
[11], but the exact value of the amplitude bi has been obtained only here. 

Why wonder about further corrections? The main reason to care about them is 
that they give rise to non-analytic terms in the momentum distribution function at 
k — K+nkp whose exact structure is characterized by the coefficients bm- Furthermore 
they give visible corrections to the off-diagonal correlators if we carefully search for 
them. This is elucidated in the there figures 2, 3, and 4 where we report p{x) for 
K = 0.1 and N = 81, for k ^ 0.5 and N = 101, k = 0.9 and N = 121 respectively. 
In the three figures, the top panels show the real (left) and imaginary (right) parts of 
p{x) directly calculated from the Toeplitz determinant. In the second panels from the 
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Figure 3. OfF-diagonal correlation function for ft = 0.5 and A'^ = 101. 
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Fourth rows: Further subtraction of the analytic contributions di (right) and 
ai,ci (left). Bottomost panels: Further subtraction of the analytic corrections d2 
(left) ai and ci (right). 



top we report real and imaginary parts of 



p{x) 



1 



p{x) [N sin Trx/L]^' 



1 , 



(51) 



Pfh{x) 

that should tend to zero for — > oo for any k ^ 1 (at k = 1, the first harmonic 
6„i is of order 0(iV")). We can in fact see that these are much smaller than the 
leading contribution for k far from 1, but become very relevant at k = 0.9 also at 
N = 121. Furthermore the second panels show very interesting oscillating behavior 
due to the first harmonic term. Thus, in order to kill the first oscillating behavior 
we subtract from pc{x) the first oscillating terms in Eq. (50), i.e. the terms in b±i. 
This subtraction is shown in the set third (from top) of panels. The effect of this 
subtraction is rather peculiar. In fact while for k = 0.1 it mainly kills the oscillations, 
leaving the absolute value almost unchanged because of the smallncss of the exponent 
in b±i, for larger k the absolute value goes down drastically, but strong oscillations 
are still present. The left oscillations are the second harmonics 6_2. In the imaginary 
part of what is left, one can clearly recognize a l/tauTrx/L shape due to the term 
di for any value of k. In fact, by subtracting the c?i term, as shown in the right 
fourth (from top) panels the value goes still drastically down of one- or two-order of 
magnitude. Oppositely the real part at the third level has a shape depending on k. 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



x/L x/L 

Figure 4. As Fig. 3 with k = 0.9 and N = 121. 

This is very easily understood: For k < 1/2 the larger term left in the expansion is 
c?2, going always like iV~^, while for k > 1/2 it is the term in ai going like 3 — 2k. At 
K = 1/2 the two terms are both N^'^, but being c?2 very close to zero (see Fig. 1), 
the most important term is ai. For these reasons, in the left fourth panel is shown for 
K = 0.1 the subtraction of ^2 and in the other two cases the subtractions of ai and ci. 
Evidently these subtractions are the correct ones and in fact decreases the considered 
value of at least one order of magnitude. Finally in the bottomost panel we report the 
final subtractions left of ai, ci or c?2 depending on the case. A further reduction of the 
value is observed. Clearly this procedure of subtraction can be repeated ad libidum, 
obtaining better and better agreement, but we preferred to stop at this level. 

A last comment concerns the absolute value of the real part of these corrections, 
that is small, but non- vanishing, as evident from the plots. The 1/iV^ contribution to 
this term can be easily obtained from second order perturbation theory, as done for 
bosons in Ref . [28] , but the asymptotic nature of this expansion makes this calculation 
not very useful. 



8. The general anyonic Calogero Sutherland model 

The Hamiltonian of the Calogero-Sutherland model is given by a long-range pair 
interactions controlled by the coupling parameter A: 

^ n2 2 / r 2 
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The ground state wave function of the Hamiltonian (52) was found in [44] and can be 
written in the form (3) with f{y) = \y\^ with normaUzation constant 

This model has always been considered only in the case of bosonic and fermionic 
statistics, or in the anyonic ones arising from having an analytic ground-state function, 
that implies k = X. While thermodynamic properties and diagonal correlations do not 
depend on k, but only on A, off-diagonal ones are a signature of the anyonic statistics 
and we can consider general k ^ A as pointed out by Girardeau [7]. The off-diagonal 
correlators are easily obtained with the same replica trick as before. 

The one-body density matrix {t = 2'7tx/L) for A'' -|- 1 particle in a ring is 

where the average is defined as 

(/ (e^^S . . . , e^^-)),_, = gii±^ ^ |A.({e'^'})|^^ / (e'^S . . . , e^^~) , 

and A]\[{{zi}) is the usual Vandermonde determinant. 

To calculate the average in (54) we use the replica trick, along the lines of the 
calculation in [27] and [28], modified to include the anyonic statistics. Namely, consider 
the function 

= (n(l - e^'0'"^(e^* - e^^O'"^^ , (55) 



N,\ 



but now the replica limit is 2np — > A. The duality transformation [34], which enables 

to re-cxpress the A^-dimensional integral (55) depending on the parameter 2np as a 
2np-dimensional integral depending on as a parameter is the same as for bosons 
[30] 

^2"pW = ^^T7IT / rf'"^^|A2„p(K})l^n[^«(i-^«)]^"'[i-(i-e'*K]^.(56) 

The duality A 1/A becomes evident by comparing the power of Vandermonde 

determinants in the previous equations. We put -^2np(0) = 1, so the normalization 
constant is given by the Selberg integral 

/■I 1 2"^' r2 (2.) Y (^ + ^] 

"'O o=l a=l ^ + A/ ^ V A / 

One can perform the sum over the saddle points at x± = 0, 1 arriving to the same 
form as before 

2np TTl it{N+2np/\)(np-l) 

^2np(t) - ij (-o vMi2 + (2np-0')A ' ^ ' 

1=0 ^ > 

but with a A dependent H\^^{X) factor given for integer 2np = m in Ref. [30] as 
TTl,y,_ r(m+l) /KA)/^-,(A) 

r(; + i)r(m-; + i) ^^(i/a) ' ^ ^ 




Figure 5. One-body density matrix gi{x) for TV = 51, for values of A and several 
values of k (that are the same in the four plots). We consider the sum of the first 
few harmonics: 5 harmonics are considered for A > 1 and only one for A < 1, 
because as-well known, the asymptotic feature of the series are more severe for 
these value of A. Distances x are measured in unit of the mean density. 



with S',„(l/A) given above and //(A) 
After trivial algebra we get 



(nLirKA]r[i + a/A])/r[i + i/A]' 



<(A) = n 



r(f) 



nr( 

a— 1 ^ 



-I 2 



A ' 



(59) 



For A = 1 this formula reduces to Eq. (25) whose analytic continuation to complex n 
can be written in terms of Barnes G-functions. Unfortunately this is not possible for 
general A (special formulas can be found for integer A or 1/A, but they are not very 
illuminating). Not having this nice representation in terms of special functions, the 
analytic continuation becomes less simple, but still straightforward and it is reported 
in the appendix. 

The anyonic statistics is obtained by setting np — I = nq — m and so 

oo rr>'/2-K/2-»n/» \ i,rx(K+2m)(7V+l)/L 
V f_1l('=/2+'»)£A 



m— — OO 



{2N\ sin7ra;/L|)(V2+2(m-HK/2)VA) 



where 



l/2-K/2-m 



(A) is intended to be continued as explained in the appendix. 
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Explictly we get 



2A/2+kV(2A) a 



Y(i J 2 I m+l \ 

V2 2A ^ 2X > 



■p/'l_ _K_ rn_\ 

2\ 2\> 



2 



22(2to+2k+1)/A j-g-|^-j 



and 5o = 1. This expression agrees with the general harmonic expansion (34) with 
K = as well known (the parameter K do not depend on the statistic k and so 
it is the same derived by Ha [45] for k = A). In Fig. 5 we report the correlations 
obtained by summing few of the first terms in the harmonic expansion for several 
values of the anyonic parameter k at fixed A to show the crossover from the fermionic 
to the bosonic curves. The extreme results k = 0, 1 were already determined in Ref. 
[30] (small manipulations are needed to show that they are identical), where a perfect 
agreement with quantum Monte Carlo simulations has been shown. 

This harmonic expansion has been previously calculated in specific cases (bosons 
and fermions for generic A in [30] and for k = A in [45]). It is a long, but simple 
calculation to show that all these expansions agrees. Other exact results appear in 
Refs. [46, 47, 48, 49, 50, 38, 51, 52] also for dynamical correlation functions. 



9. Conclusions 



In this manuscript we proposed a variation of the replica method that allows the 
calculation of the large distance asymptotic of off-diagonal correlation functions in 
anyonic models. The main result we found is that the replicated correlation Z2np{t) 
is a function with 2p branches corresponding to different anyon statistics k = q/p 
with q = 0, . . . 2p — 1. After having understood this, it is almost straightforward but 
tedious, to obtain all the harmonic terms of the large-distance expansion from known 
results for bosons. We applied the method to a gas of impenetrable anyons (Tonks- 
Girardeau) and the most general anyonic Calogero-Sutherland. For the former a 
number of analytic and numerical checks confirmed the validity of the replica approach. 
For the latter instead, for special values of the interaction parameters (i.e. integer or 
rational), a number of exact results have been already found thanks to the random 
matrix theory and its connection with Jack polynomials [47, 48, 49, 45, 50, 38]. In some 
cases the correspondence of the two approaches is obvious, while in some others (for 
rational \ = p/q the correlation functions are given by sum over fractional excitations 
involving p + g quantum numbers, giving an irregular function of A with no clear limit 
to irrational values) it impHes particular relations for Jack polynomials. Showing in 
full generality this correspondence remains an open problem. 

At the same time, another open problem is the correspondence of our approach to 
the one by Patu, Korepin and Averin [22] based on Fredholm determinants. We stress 
that (as far as we are aware) the explicit correspondence of the two approaches has 
not yet been proved even for bosons. The results coming for these two sets of works 
give an almost complete physical picture of the behavior of the off-diagonal correlator 
of the gas of impenetrable anyons. Beyond the impenetrable limit, only the Bethe 
solution of the Lieb-Liniger model is known [4, 5], with still no attempt to tackle the 
hard problem of the correlation functions exactly. Only the power-law structure (and 
the corresponding singularities) are known from bosonization [11, 12]. A powerful 
approach would be to generalize the quantum inverse scattering methods for bosons 
[53] to anyons and then to mix integrability and numerics to get the full correlation 
functions from the form- factors (on the line of Ref. [54] for bosons). 
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Note Added 

After the completion of this work, a manuscript appeared where some large distance 
asymptotic of the off-diagonal correlations for the Tonks-Girardeau gas have been 
found at finite temperature [55]. The results obtained there arc different and 
complementary to ours. Notice that the integral system of partial differential equations 
they derive is independent on the anyonic parameter k that enters only through the 
initial condition, on the line of our previous [21] and current results. 
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Appendix A. The analytic continuation 

We need to analytic continue terms of the form 

r(a/A) 



hi 



n 



lJ.r{{m-l + a)/X)- 
Using the integral representation of the log of the Gamma function 



iogr(z) = 



dt 

T 



+ {z- l)e- 



(A.1) 



(A.2) 



1 - e-* 

summing the resulting geometrical scries and changing integration variable, we have 



log hi 



Jo t 



1 - e-" - e-*(™-') + e-"* (m - 1)1 ^-(a-i)*' 



(l-e-*)(l-e-^*) 



A 



(A.3) 



For the first harmonic and any k and A, the single terms in the integral all diverge, 
only their sum is finite. The integral is an analytic function of I and m and can then 
be used to evaluate /i^ for any complex to, I for the first harmonic. The other are 
obtained from the recurrence relation in the text. 

For the coefficients iJ^(A) in the text we also need the analytic continuation of 
fm = r((TO + c)/A)/r(c/A) that are very similar 



log/w 



f 

Jo 



dt 



(1 



-mt\2 



m 



-(A-l)f 



(A.4) 



.(l-e-*)(l-e-^*) A 
that is vanishing for A = 1 and m — > 1 as it should. Plugging everything together 

log Hl{\) = 2 log hi + log fm . (A.5) 
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